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1.0 SUMMARY 


This paper presents the results of a study to determine the source of 
computational error in the IMU onorbit alignment software as specified by the 
Level C IMU SOP FSSR. Simulation runs were made on the IBM 360/70 computer 
with the IMU onorbit alignment software coded in HAL/S. The results indicate 
that for small IMU misalignment angles (less than 600 arc seconds), single 
precision computations in combination with the arc cosine method of eigen 
rotation angle extraction introduces an additional misalignment error of up 
to 230 arc seconds per axis. Use of the arc sine method, however, produced 
negligible misalignment error. As a result of this study, the arc sine 
method has been recommended by means of a software change request (CR) for 
use in the IMU onorbit alignment software. This CR has been approved for 
STS-1 and beyond. 


2.0 INTRODUCTION 


During IBM's Level 2 IMU onorbit alignment testing, it was first suspected 
that the largest source of software produced alignment error for small 
torquing angle misalignments might be the flight control utility routine 
MAT_TO_QUAT. After review of these initial test results, it was decided 
by JSC, IBM and MDTSCO that detailed studies of this problem be made at 
MDTSCO to confirm the initial conclusion. The proposed studies were aimed 
at determining whether the additional computational error was in fact 
caused by MAT_TO_QUAT or by the IMU alignment software design. The questionable 
design centered around the use of single precision computations in combination 
with the arc cosine method for extracting the eigen rotation angle from 
the error quaternion, from which torquing angles were computed. The results 
from these studies indicated that tl... iatter was the cause of the error. 

Based on this conclusion, change reques* IS123A '•'as submitted and approved 
to change the eigenaxis rotation am,, le extraction method from the arc 
cosine method to the arc sine method for STS missions number 1 and beyond. 

Section 3.0 contains a discussion of the equations which define the IMU 
onorbit alignment software and a description of the simulation runs actually 
made. Section 4.0 contains a discussion of the results, and Section 5.0 
presents the conclusions and recommendations. 
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3.0 DISCUSSION 


The IBM 360/70 computer was used in simulating the HAL/S version of the 
IMU onorbit alignment torquing angle software for this study. This computer 
was used for the following reasons: (1) the IBM 360/70 computer implements 
a HAL/S compiler which enables the simulation software to be coded in 
UAL/S, ( 2 ) this computer emulates the word size of the AP101 flight computer. 
Although the word size is equivalent for both computers, AP101 double 
precision computations are not as accurate as the IBM 360/70 computer; 
therefore, the double precision results presented in this study are not 
representative of AP101 double precision accuracy. 

The following sections describe the simulation software and test runs 
for determining the source of computational error in the onorbit alignment 
software design. 


3.1 DETERMINATION OF THE MEASURED PRESENT-CLUSTER-T0-M50-MATRIX 


In order to determine the measured transformation matrix from present 
cluster to M50 coordinates, the onorbit alignment software obtains measurements 
of two star LOS (Line of Sight) unit vectors in present cluster coordinates 
by means of the STAR_TRACK or the COAS SIGHT mode, that is 


^m ” ^xm 1 pc 
^m = ^xm^pc 


+ s ymJpc + s zm^pc (1) 

+ ^ymjpc + ^zm^pc (2) 


In addition, the computer has stored in memory the exact LOS unit vectors 
in M50 coordinates of these two stars as follows: 


\ = Sxe^SO + s yeJ'M50 + s ze k M50 (3) 

T e = T xe i M50 + T ye3M50 + T ze%0 ( 4 ) 


These two sets of unit vectors, measured and exact, are used to compute 
the desired transformation matrix. 

First, the exact LOS unit vectors given in M50 coordinates are used to 
form three unit vectors defining a star coordinate system expressed in 
the M50 coordinate system as follows: 


Uxe 

- s e 

(5) 

Uye 

= UNIT (S e x T e ) 

(6) 

*^ze 

= UNIT (U xe x U ye ) 

(7) 
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Next, the three unit vectors are used to form the columns of the transformation 
matrix from star to M50 coordinates: 



u xc (l) 

Uye {1) 

U ze (l) 

U xe (2) 

U ye (2) 

U ze (2) 

_U xe (3) 

U ye (3) 

U ze (3_) 


( 8 ) 


In a similar manner, the measured LOS unit vectors (given in present cluster 
coordinates) are used to form this same star coordinate system expressed 
in the preset. t cluster coordinate system as follows: 


i = S 
xm J m 

(9) 

iym = UNIT (S m x T m ) 

(10) 

'ztr ~ UN 11" (U xm x U ym ) 

'll) 


These unit vectors are then used to form the rows of the measured* transformation 
matrix from present cluster to star coordinates: 


Hr "I 


Uxm( 

U xm (2) 

U X m(3) 

L t pcJ 

= 

Uy m ( 1 ) 

Uy m (2) 

Uym( 3 ) 



_Uzm( 1 ) 

U zm (2) 

Uzm(3)_ 


Finally, the measured transformation matrix from present cluster to M50 coordinates 
is formed by the following matrix multiplication: 


^M50 
T PC = 


M50l|^S ; 

Js JL t pcJ 


(13) 


3.2 COMPUTATION OF THE MISALIGNMENT MATRIX AND ASSOCIATED ERROR QUATERNION 

The misalignment matrix is computed as a function of the stored transformation 
matrix from M50 coordinates to desired cluster coordinates (desired REFSMMAT) 
and the measured transformation matrix from present cluster to M60 coordinates, 
that is 


'vDC 

T PC 


~ DC 1 

kM50l 

_ t M5oJ 

! JPC J 


(14) 


* The tilde over the matrix indicates that the matrix is derived using onboara 
measured quantities and is to be distinguished from matrices without the 
tilde which are considered exact. 
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The error quaternion AQ is defined by the followup convention 


~AQS* 


cos iti/Z 

_AQV_ 


u sin uj/2_ 


(15) 


where AQS is the scalar part of the quaternion, AQV the vector part, u> 

A 

the eigen rotation angle about the eigenaxis, and u the eigenaxis unit 
vector. In order to compute the error quaternion associated with the 
misalignment matrix (14), the utility routine MAT_TO__QUAT is called, that is 


CALL MAT_TO_QUAT (T pc ) 

ASSIGN (AQ) 

This error quaternion is the basis for determining the torquing angles 
for platform realignment. The eigenaxis direction is obtained by extracting 
a unit vector* from the vector part of AQ. The eigen rotation angle can 
theoretically be obtained either from the scalar part of AQ or the magnitude 
of the vector part of AQ. The methods of extracting the eigenaxis rotation 
angle are described in more detail in the next section. At the time this 
study was performed, the eigen rotation angle was extracted by use of 
the scalar of the error quaternion. 


3.3 DESCRIPTION OF EIGEN ROTATION ANGLE EXTRACTION METHODS 

There exist two methods for extracting the eigenaxis rotation angle: 

(1) Arc cosine method (using the error quaternion scalar) and (2) Arc 
sine method (using the error quaternion vector). The equations for both 
methods are now presented. 


3.3.1 Arc Cosine Method 

The first method initially used by the IMU alignment software design is 
the arc cosine method. This method uses the scalar part of the error 
quaternion AQS from equation (15), that is 


cos w/2 = AQS 


(16) 


* Although this paper does not discuss problems associated with this extraction 
process, it should be mentioned that for values of misalignment eigen 
rotation angles less than 10 arc seconds, the unit vector associated 
with the eigenaxis direction can be in error as much as 100 degrees. 

However, for such small eigenaxis rotation angles the direction of rotation 
is irrelevant. 


After taking arc cosine of both sides of this equation and multiplying 
by two, the eigenaxis rotation angle is given by 

a) = 2 cos (AQS) (17) 


3.3.2 Arc Sine Method 

The second method for extracting the eigenaxis rotation angle is the 
arc sine method. This method uses the magnitude of the vector part of 
the quaternion AQV from equation (15), that is 


s i n a>/2 = l*Qv! 


(18) 


After taking arc sine of both sides of the equation and multiplying by 
two, the eiproaxis rotation angle is given by 


w = 2 sin" 1 |2QV| 


(19) 


3.4 DESCRIPTION OF TEST SETUP AND RUNS 

For simulation purposes, the exact LOS unit vectors are arbitrarily defined 
as follows: S e is defined by an azimuth angle of 60 degrees and an elevation 
angle of 30 degrees, Tp is defined by an azimuth angle of -60 degrees and 
an elevation angle of 30 degrees. The measured LOS unit vectors are derived 
by rotating the exact LOS unit vectors (S e , f e ) from the M50 frame to 
actual platform coordinates, that is 


DC 

M5C (MEASURED) 
_M50( ACTUAL) 

^ i 

M50 ^ 

5 e 

DC 

~M50 (MEASURED*) 

T 

M5oJ 

_M50 (ACTUAL) 

‘e 


( 20 ) 


( 21 ) 


The input misalignment eigenaxis direction is defined by an azimuth angle 
of 45 degrees, and an elevation angle of 30 degrees. The eigenaxis rotation 
or input misalignment for each specific test case is given in Table 1. 

For each case, an input misalignment quaternion is formed from these two 
input quantities and then a call to QUAT_TO_MAT results in the input misalign- 
ment matrix. Also, for each case, the measured and exact LOS unit vectors are 
unitized to IBM 360/70 single and double precision since they were generated 
using HP9825A precision. The stored REFSMMAT is defined by a quaternion which 
has an eigenaxis direction and eigenaxis rotation angle. The eigenaxis direction 
is arbitrarily defined by an azimuth angle of 45 degrees and an elevation angle 
of -3U degrees. The eigenaxis rotation angle between the M50 frame and desired 
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X, 


cluster coordinates is taken to be 90 degrees. A call to QUAT_T0_MAT results 
in the stored transformation matrix from M50 to desired cluster coordinates. 

For each case, this REFSMMAT is orthogonal i zed in IBM 360/70 double precision 
before being used in th« onorbit alignment flight software. 

Each test run used the onorbit alignment software design as described 
in sections 3.1 and 3.2. For each input misalignment angle case, both 
single and double precision versions of the HAL/S code were exercised 
(Appendix A and Appendix B, respectively) and an error computed by taking 
the difference between the input misalignment angle and the output misalignment 
angle for both eigenaxis rotation angle extraction methods. This error 
became the standard of measure for determining the best method. 
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4.0 RESULTS 


A summary of the IMU onorbit alignment HAL/S simulation results is presented 
in Table 1. The output misalignment angle and the misalignment error 
Sid are given as functions of the input misalignment angle <j, for both 
precision versions (single, double) and both angle extraction methods 
(arc cosine, arc sine). The range of coj is from 0 to 648000 arc seconds 
(180 degrees). As inoicated in the last column of the table, both extraction 
methods result in zero misalignment software error over the entire range 
of input rotations when double precision is used. The single precision 
error column of &o indicates that the arc cosine method has a sizeable 
error for small misalignment angles. In order to visually show the amount 
of misalignment error for both methods when codec in single precision, 
a graph of misalignment error versus input misalignment urj for the 
arc cosine and arc sine rotation angle extraction methods is presented in 
Figure 1. For the arc cosine method, Figure 1 shows approximately 400 
arc seconds of misalignment error for input misalignment angles near 0 
arc seconds. This 400 arc seconds of misalignment error would result 
in approximately 230 arc seconds per axis of additional misalignment error 
to the current IMU alignment error budget. For increasing input misalignment 
angles, a gradual decrease in misalignment error can be seen. For the 
arc sine method. Figure 1 shows that for input misalignment angles near 
C arc seconds, the misalignment error is negligible (i.e., .03 arc seconds). 
For large input misalignment angles (near 648000 arc seconds), the maximum 
misalignment error is 340 arc seconds. This amount of misalignment error 
would result in approximately 196 arc seconds (per axis) of additional 
misalignment error to the IMU alignment error budget. 
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5.0 CONCLUSIONS AND RECOMMENDATIONS 


For small IMU misalignment angles (less than 600 arc seconds), single 
precision computations in combination with the arc cosine method introduces 
additional misalignment error of up to 230 arc seconds per axis. The 
arc sine method, on the other hand, produces results equivalent to double 
precision accuracy. The misalignment error 6u>, obtained during this study, 
is clearly a function of the eigenaxis rotation angle extraction method 
and the precision used in the HAL/S code. The presence of MAT_TO_QUAT 
did not significantly influence the amount of misalignment error. 

It is therefore recommended that the arc sine method of eigen rotation 
angle extraction be implemented into the onorbit alignment software at 
the earliest possible date. This software change would eliminate a known 
computationally-induced bias error. 












TABLE 1.0 - SUMMARY OF I MU 0N0RB11 ALIGNMENT HAL/S SIMULATION RESUL1S 


INPUT MISALIGN 
MENT ANGLE, Wf 


(arc seconds) 

0 


OUTPUT MISALIGNMENT ANGLE, <o 0 
_ ( arc seconds) _ ^ 

8 Engle Precision - ' Double' Precision 


MISALIGNMENT CRROR, 6u 

_(arc _s> 

S i nql e Prec ) sTon 


(arc sec onds) _ 

Double Precision 


15" 
" 25" 

50 
75" 
TOO 
1 75 
150 

m 

300 
600 
900 
1800 
2700 
“ 1600 
640800 
644400 
646300 
646200 
647100 
647400 
647700 


co's^ 1 * 
si n _1 * 
cos » 
si n' 1 * 
cos * 
s i n ' 1 - 
cos 
si n 
cos 
sin 
cos 
s i n 
cos -1 
sin 
cos 
s I n 
cos 
sin 
cos 
s i n 
cos 
s i n 
cos 
s i n 

COS"' * 

s i n 
cos 
s I n 

cos “4 - 
sin 
cos 
s i n 
cos' 1 * 
sin -1 * 

-U 


1, 
-1. 
1: 
-1, 
-1- 
■1, 
•1 

•1- 
; - I 

■ u 

■1. 

■ 1 - 
i" ' - 

: - T » 
•1. 

1 

■1 

•-U 

;-u 

•1 

-u 
■ 1, 
• 1- 


COS" 1 * 


s i n 
cos 
s i n 
cos 
s I n 
cos 
s i n 
cos 
sin 
cos 
sin 


-l s 
-L 
- T- 
- 1 . 
-1. 
;-i 
i- 
i. 
-i- = 
i. 


407.86 

'o.o " 

COS'** 

mar 

070 

0.03 

0.0 

sin" 1 * 

0.03 

0.0 

40? .86 

5.0 

cos' 1 * 

397.86 

0.0 

5.00 

5.0 

sin' 1 * 

0.00 

0.0 

407". 85 

" 15. XT' ' 

cos' 1 * 

137 ITT 

0.0 

15.01 

15.0 

sin -1 * 

0.01 

0.0 

407.86' 

. - 

COS'] * 

"377186“ 

0.0 

?5.02 

25.0 

sin' 1 * 

0.02 

0.0 

407.86 

50.0 

COS'" 1 "* 

157:86' 

~~o 

50.0? 

50.0 

$ i n “ * * 

0.02 

0.0 

407.86 

75.0 

COS' 1 = 

377.86 ' 

1.0 

74.99 

75.0 

sin -1 = 

0.01 

0.0 

403.86 

100.0 

COS' 1 * 

302.86 

O'.f) ' 

99.99 

100.0 

sin' 1 * 

0.01 

0.0 

407.86 

' 125.0 

cos' 1 •- 

277.86' 

"0.1" " ' 

175.01 

125.0 

sin' 1 * 

0.01 

0.0 

407.86 

150.0 

COS' 1 - 

752 .86 

0.0"" ' 

150.0? 

150.0 

sin' 1 * 

0.0? 

0.0 

407.86 

200.0“ ' 

. 

cos" * * 

202.86“ 

0.0 

700.00 

200.0 

sin' 1 * 

0.00 

0.0 

569.73 

300.0 

COS' 1 * 

7655773 " 

- " 0:1 

300.07 

300.0 

s i n' 1 = 

0.02 

0.0 

697./ 7 

600.0 

COS' 1 * 

‘97.77 ' 

0.0 

600.01 

600.0 

sin' 1 * 

0.01 

0.0 

986.80 

900.0 

COS' 1 = 

86.80 

0 - 

900.00 

900.0 

sin' 1 * 

0.00 

0.0 

1846.14 

1800.0 

COS' * 

4(7.14 

179" 

1800.01 

1800.0 

sin' 1 = 

0.01 

0.0 

273?'. 33 

" 2700.0 

COS' 1 - 

32733 

0.0 

2/00.01 

7/00.0 

sin' 1 * 

0.01 

0.0 

"' 3625.75 

3600.0 

COS' -1 * 

"75'. 75 

0.0 

3600.01 

3600.0 

sin' 1 * 

0.01 

0.0 

640799.75 

640800.0 

COS' * * 

0.25 

1.0' ' 

640/90.43 

640800.0 

sin' 1 * 

9.57 

0.0 

644399.56 

644400.0 

COS' 1 * 

0.44 

0.0 

6443/6.81 

644400.0 

sin' 1 - 

23.19 

0.0 

645299. 75 

645300.0 

cos' 1 = 

0.25 

o;o 

64576/. 50 

645300.0 

sin- 

32 . 50 

0.0 

646199.87 

646200.0 

cos' 1 - 

0.13 

1.1 ' ' ' ' 

646153.50 

646700.0 

sin' 1 * 

46.50 

0.0 

647099.62 

647100.0 ' ' 

cos' * 

0.38 

0.0 

64/073.31 

64/100.0 

r in' * 

76.69 

0.0 

647399.75 

647400.0 

cos' 1 * 

0.25 

"O.o 

64/246.06 

647400.0 

sin -1 * 

153.94 

0.0 

647699.68 

647700.0 " 

COS' 1 = 

0.32 

0.0 " 

647527.31 

647700.0 

sin' 1 * 

172.69 

0.0 
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TABLE 1.0 - SUMMARY OF IMU ONORBIT ALIGNMENT HAL/S SIMULATION RESULTS - Concluded 


INPUT MISALIGN 
MENT ANGLE, 10 , 
(arc seconds) 1 

- OUTPUT MISALIGNMENT ANGLE, u> 0 
(arc seconds) 

MISALIGNMENT ERROR, 
(arc seconds] 

, 6u) 

! „ , 

Single 

Prec ision 

Double Precision 

5ingle 

Precision 

Doub 

le Precisii 

647800 

cos‘|= 

sin'J= 

" 6477 ' 98 . 8 'r 

647549.37 

647800.0 

647800.0 

cos'|= 

sin*|= 

0707 

250.63 


U 70 

0.0 

647850 

COS‘|= 

sin”J= 

647848.81 

647549.37 

647850.0 

647850.0 

cos - |= 

sin'f= 

0718 

300.63 


UTO 

0.0 

647870 

cos" |= 

si n" 

647869.56 

647549.37 

647870.0 

647870.0 

cos‘j= 

sin’ 1= 

0.44 

320.63 


0.0 

0.0 

”647880 

cos']= 

si n“4 = 

647549 !37 

647890. 0 
647890.0 

COS"J= 

sin _ l= 

0.32 

340.63 


0.0 

0.0 

647815 

COS _ l= 

sin"l= 

647814.75 

647715.00 

647915.0 

647915.0 

cos*|= 

sin'*= 

0.25 

200.00 


0.0 

0.0 

647965 

cos*l= 

sin -1 = 

647964.62 

647715.00 

647965.0 

647965.0 

COS‘j= 

sin"l= 

0.38 

250.00 


0.0 

0.0 

648000 

COS‘|= 

sin -1 = 

647898.81 

647753.18 

648000. 0 
648000.0 

COS'J= 

sin _ l= 

0.19 

246.82 


0.0 

0.0 
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APPENDIX A 

(HAL/S CODE IN SINGLE PRECISION) 
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